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ABSTRACT 

The full treatment, both nonperturbative - through 0{l/m? c ) - and perturbative - through 
0(a s ) and BLM type 0{a 2 8 ) - of the semileptonic lifetime of the D meson is given. A dedicated 
discussion of the numerics involved is made, with the result that the leading order prediction 
plus corrections fail to match experiment. To explain the difference, a model of duality is 
invoked, where it is shown that the asymptotic nature of the nonperturbative expansion may be 
responsible for the theoretical experimental discrepancy. A particularly interesting conclusion 
which follows from the model is that any attempts to extract a s from t decays are doomed, 
since the same duality violating mechanism is at work here. 



i 



Introduction 



"All is mystery; but he is a slave who will not struggle to penetrate the dark veil" - 

Disraeli 

Hopefully, there will be a few people outside the high energy theory communitity 
who will at least take a glance at this thesis. This introduction is written for them, 
in hope that all the complicated equations, diagrams, and verbage following will have 
some general meaning. I think it was Feynman who said that you should be able to 
explain your research to a seventh grader, or something like that. How to DO the 
research, is, of course, another matter all together, but what the problem is, and what 
its solution is should be understandable to anyone. This is my attempt to explain: 

Most people know that there are things called atoms, and have seen simple diagrams 
(resembling the solar system to some extent) which picture electrons whirling around 
a nucleus. Although, after a semester of quantum theory, one would most likely cringe 
at such a diagram, these simplified pictures do in fact illustrate some key features 
of the atom. The electrons are usually displayed either as clouds, illustrating the 
wave nature of the very small, or as particles traversing orbits around a center, the 
nucleus. Usually, the electron, when pictured in the 'planet' like orbital picture, is 
shown as a 'dot'. In some sense, this is a fair representation. As far as we know at 
present, the electron is a particle which has no structure - no radius whatsoever, and no 
constitutents. The electron is, as far as we can tell, the perfect physical representation 
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of a mathematical point. The nucleus, however, is an entirely different story. It has 
constituents, being made up of protons and neutrons. Surprisingly, although protons 
are positive in charge (neutrons are neutral!) the nucleus does not 'fly' apart (like 
charges repel). There is another force binding the nucleus together, stronger than the 
electrostatic repulsion, called the strong force. The modern day theory of the strong 
force is called QCD - Quantum Chromodynamics, which is the topic of my thesis. 
The 'objects' of QCD, however, are not the proton and neutron. It is known for some 
decades now that the proton and neutron are not at all like the electron. These particles 
do have internal structure. Indeed, the proton and neutron (and many other particles 
like them called either mesons or baryons) are made from particles called quarks and 
gluons - the true 'objects' of QCD. Quarks and gluons themselves are like the electron 
- seemingly pointlike. Quarks come in different 'flavors' (i.e. there are different kinds 
which differ by their masses). So far, there are six different flavors found: up, down, 
strange, charm, bottom, and top. Quarks are never found alone, but are confined either 
in pairs (mesons) or triplets (baryons). As stated above, there is a theory - QCD - which 
mathematically details the dynamics of how quarks interact. What does it mean to say 
that there is a 'theory' detailing the dynamics? It basically means we should be able 
to calculate some property of quarks, and compare it to experiment. For example, this 
thesis is concerned with calculating the lifetime (the strange, charmed, bottom, and 
top quark are all unstable, and quickly decay when created at a particle accelerator) of 
a a meson called the D meson (the D meson is a pair of quarks - one charmed and one 
up or down). Specifically, I calculate the decay of the D meson into a pair of leptons 
(an electron is a member of the lepton 'family') and all possible hadrons (a hadron is 
either a meson or a baryon, but in general I mean here all other possible combinations 
of quarks). The D meson lifetime is well established experimentally. Calculating it, on 
the other hand, is absolutely difficult. QCD is a notoriously difficult theory to calculate 
with. Noble Prize winner Steven Weinberg is quoted as saying: 
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"There's a long tradition in theoretical physics, which by no means affected everyone 
but certainly affected me, that said the strong interactions are too complicated for the 

human mind." 



Indeed, making an exact calculation of the decay time of the D meson, or any other 
meson will probably never happen. This is really not such a big deal, however, as this is 
almost always the case in physics anyway. Hardly anything can be calculated without 
making some approximations, and when it can, it usually doesn't correspond to the true 
physical world anyway. (Of course, nature doesn't depend on whether we can 'calculate' 
it or not - e.g. atoms do not calculate their own behavior!, still much of the richness of 
our universe resides in the fact that nature is calculable only in approximation. If nature 
happened to correspond to a very simple theory which could be used to calculate for all 
physical phenomenon, we would live in a boring place indeed. (Let me remark that it 
seems certain that nature abides by some very elegant, and simple theories which lead 
to, most of the time, very rich phenomena. It is usually very 'easy' to write the theory, 
and most often very difficult to calculate with it - thus the need for approximations!)) 
The D meson is seemingly an ideal situation for making approximations. The reason is 
that the charmed quark has a high mass. What can be done then computationally, is 
to express the decay lifetime in a power series: 

r = a + + a 2 (^) 2 + a 3 (^) 3 ... (1) 

m c m c m c 

here m c is the mass of the charm quark, and fi is a tool of the theorist, and is chosen so 
that the expansion converges (there are subtleties here which even real experts do not 
get!). Numerically, fi is chosen so that fj,/m c = .5 or so. Suppose then, that ao, ai, a-2, ■■■ 
are - all different - but all close to each other (unless there is some strong reason, this 
could be expected - a statement of vague intuition to be sure). Then, for each new term 
in the series, say 04, 05, etc. the new term is also multiplied by increasing powers of 
fi/m c = .5. Technically, it is impossible to calculate the full series - no one knows how, 
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but we can approximate the result by calculating the first few terms and then truncating 
the series when each new term gets insignificant (or too hard!). The calculation of ao, 
a±, and 02 were all completed before I knew what a heavy quark even was, but not 
much before! Check out the result of this calculation - equation 57, and in particular 
note that m c is in the denominator of the last few terms, in tune with what is written 
above. Essentially, the first five chapters (aside from the preface) are concerned with 
this computation. All the details of 'how to do it' are here. Chapter 6 and 7 explain 
the actual numerical values that are used to check and see if equation 5.22 actually 
matches experiment. This is an interesting area of scientific debate. As it turns out, 
if one measures numerical values for the charmed quark mass and other parameters 
from experiments and applies them to the problem of the D lifetime, the theoretical 
value of the decay time is roughly 50 percent of the experimental value (pretty bad 
agreement) - check out equation 7.3. The key numerical parameter here is m c , since, 
as one can see from equation 5.22, it enters into the expression of the lifetime as m^ c . 
Because m c is to the fifth power, only a slight increase, roughly 10 percent, of it's 
value will compensate the 50 percent deficit between theory and experiment. However, 
those (Dikeman included) who think that m c should have a low value and thus spoil the 
agreement between theory and experiment regarding the D lifetime believe this because 
it would spoil too many other results which say m c is low. Thus, there is a problem. 

One possible next move is to try and calculate the next term in the series and see 
if it somehow can raise the decay value a whole 50 percent - an unlikely scenario since 
remember, each new term should be .5 times smaller than its previous one, and thus 
numerically - check out equation 88 - we would expect, at best, this 'cubic' term in the 
series to contribute at most 10 percent, but not 50 percent. Still, though, one must 
try, and so the calculation of the next term in the series, 03, was my first task as a 
graduate student. It was hoped that this 'cubic' term would explain the puzzle. This is 
the discussion of chapter 8. Well, after a lot of work, no such luck. My advisor Misha 
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Shifman, his past student Boris Blok, and myself wrote a paper demonstrating that the 
next term in the series doesn't help matters, and may even make them worse. At least 
in this paper we tried to give a number of possible explanations why the puzzle existed 
- our main guess was that the charmed quark was not quite heavy enough for the power 
series expansion to work in the first place. Despite the negative results, though, the 
paper was mildly successful - currently it has been cited in other papers 17 times, which 
is not at all bad. 

With the puzzle deepened, I went on to another project, but then returned to the 
D meson puzzle a year later. My advisor, Misha Shifman, suggested to me and fel- 
low graduate student and friend Boris Chibisov while in a summer school at Boulder, 
Colorado that there may be a way of estimating what is 'leftover' when truncating the 
series, i.e. estimating the remaining pieces, a?, ... aio, etc. of the series. This work was 
the first real attempt at such an estimate concerning the D meson, or any other QCD 
phenomenon, and so is a much more general paper than the first one, and a bit more 
successful - it has 21 citations so far, and I expect a lot more in the future. The ideas 
of this paper are presented in chapter 10. In general, these ideas are interesting, but 
they don't exactly solve the remaining puzzle, instead putting it in proper perspective. 
As I allude to above, it is difficult to solve a hard problem without making approxima- 
tions, and chapter 10 is loaded with them, but we at least hope to correctly capture 
some important features of the real, if there ever really is one!, solution to the puzzle. 
In addition, the model developed in chapter 10 has strong implications for other QCD 
processes as well, and sheds a lot of doubt on a few particular phenomenological studies 
which are quite important. And thats the end of the story. 
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Chapter 1 



Preface 

"The age in which we live is the age in which we are discovering the fundamental laws 
of nature, and that day will never come again." - Feynman 

"The more the universe seems comprehensible, the more it also seems pointless. But if 
there is no solace in the fruits of our research, there is at least some consolation in the 

research itself. Men and women are not content to comfort themselves with tales of 
gods and giants, or to confine their thoughts to the daily affairs of life; they also build 
telescopes and satellites and accelerators, and sit at their desks for endless hours 
working out the meaning of the data they gather. The effort to understand the 
universe is one of the very few things that lifts human life a little above the level of 
farce, and gives it some of the grace of tragedy" - Steven Weinberg 

In any physical theory, it is often helpful, sometimes crucial and perhaps even manda- 
tory, to identify the simplest problem the theory can solve. Fermi, referring to this 
feature of theoretical physics, liked to pose the question 'What is the hydrogen atom of 
the theory?' For Quantum Chromodynamics (QCD), the theory of nuclear force, this 
is not such an easy question to answer. Historically, as a shockingly rich spectra of 
hadrons was discovered by experimentalists in the 50's, theorists attacked problems of 
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strong resonance physics not with powerful analytic results from first principles QCD - 
QCD did not even exist until 1973, but instead with ideas based on symmetries. The 
group (flavor) structure of QCD was revealed by the spectroscopy of hadrons. This is 
the work most associated with Gell-Mann and Ne'eman Gell-Mann realized that 
the spectra of hadrons - protons, neutrons, and the myriad resonances being discov- 
ered - could be explained by introducing Q what he deemed fictitious entities: quarks, 
which carried a color charge, and some 'glue' to hold the quarks together. Aside from 
Gell-Mann's 'Eight-fold way' which, anyway, dealt with the flavor chiral symmetry [|| 
of QCD, and not its more fundamental color symmetry, the answer to Fermi's question 
might be, at least regarding the gauge structure of the theory, the special puzzle of the 
A ++ wave function. 

The A ++ is a spin 3/2 particle, a fermion, but with entirely symmetric wave- function. 
As is well known, a half-integer spin particle necessarily has an anti-symmetric wave 
function, and so there is a puzzle. Instead of appealing to some perverse solution, the 
elegant notion of a color charge, yielding an antisymmetric, colorless wave function was 
put forth which solved this and other such problems [||. Thus the idea that quarks 
carried color, and appeared in nature in only confined combinations allowing a color- 
singlet wavefunction was born, and so was QCD ||, ||. In this sense, i.e. elucidating 
that there is a color wave function, and thus an SU(3) non-Abelian gauge theory, [|| 
the A ++ is the hydrogen atom of QCD. But when Fermi asked this question, surely 
he was referring to a more dynamical problem in the theory - the A ++ has helped us 
identify the theory, but now we should calculate the spectrum of the simplest system 
with it. 

Possibly the most genuinally interesting thing about QCD is that there are such 
a variety of dynamical situations where we in fact have a simple hydrogen atom like 
problem, which quickly becomes contaminated with the intricacies of QCD. Analyzing 
phenomena such as e + e~ annhilation, DIS, heavy quark decays, etc. we, at first, find 
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quick successes with traditional methods only to later suffer severe impedances where 
our methods aren't applicable. At the core of the simplicities of these systems is asymp- 
totic freedom ||. Asymptotic freedom tells us that we can use our powerful arsenals of 
perturbative methods, developed already in the study of QED, to attack problems in 
the ultraviolet (high momentum domain) where the QCD coupling constant is small. 
In the infrared (soft momentum domain), however, where a s ~ 1 there is what is called 
infrared slavery, and we must resort to other methods. Essentially, asymptotic freedom 
creates a calculational boundary, beyond which we must be extremely clever to extract 
model independant results. Due to the nature of the running of the strong coupling 
constant, we face a situation where for soft processes, typically at scales around a few 
GeV, we have no reliable way to calculate in QCD. QCD is a nonlinear system with no 
analytic solution, and in the infrared, we have no small quantity that we can expand 
in - a truly great challenge. 

And so, as it turns out, the biggest challenge of QCD is to understand its nonper- 
turbative physics. Really, then, if we are to ask Fermi's question of QCD, cast in the 
pall of trying to understand the genuine solution to the soft problem, we realize that 
perhaps we are stuck! Indeed, taking a fairly harsh view of some of the extremely hard, 
and amazing work in QCD phenomenology done in the last 20 years on the soft prob- 
lem, e.g. ref. [Jjj, one might say that the problem of the soft, nonperturbative aspects 
of QCD have, in some sense, merely been swept under the rug! Since the inception of 
QCD, the primary theoretical tool used to deal with the soft physics of QCD has been 
Wilson's Operator Product Expansion f8|. Essentially, the central utility of the OPE is 
that it allows us to characterize - but not calculate! - the fundamental soft quantities 
of QCD, which we can then extract from phenomenology. So we give a name for our 
ignorances, and can even ask for them to be measured by experimentalists, but we have 
no idea how to calculate these ferocious beasts. With that perspective, perhaps the 
answer to Fermi's question will not actually be how to calculate this-or-that process, 
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but instead how to calculate this-or-that condensate! Who knows. 

This thesis will not at all be an attempt to answer the question of how to calculate 
the condensates in QCD. I, instead, will back down from this Herculean task, and take 
on a lesser trial. In this thesis, I will look at the lifetime of the D meson. In all 
hadronic processes there is, perhaps, no system where one can more easily immerse 
oneself into the intricasies of QCD mentioned above. The D meson, in fact, is much 
like the hydrogen atom! This meson contains a heavy quark at its center with a light 
quark bound to it by an infinite number of quark/gluonic degrees of freedom. Because 
the unstable charm quark at the meson's core is so heavy, energies involved in its decay 
will be large enough to invoke asymptotic freedom, and thus traditional perturbative 
methods. However, since the charmed quark decays in an environment where typical 
momenta are of the order of Aqcd (the momenta where the strong coupling constant, 
a s ~ 1), we will also have to deal with nonperturbative effects. By solving the problem 
of the lifetime of the D meson, hopefully revealing some of nature's trickeries in the 
strong interaction section, we also may receive other dividends as well - this is a weak 
interaction, and so we might get some information on CP violation if we can do the job 
right. 

Not to spoil the story, but as it turns out, we won't be able to do the job right. No, 
in fact, we will learn from this case that the answer to Fermi's question posed above: 
'how do we calculate this-or-that condensate' still might get us nowhere. Indeed, to 
really get to the bottom of QCD we will have to know how to calculate the vacuum 
fluctuations themselves, and how this could ultimately be done, no one knows. Well, 
actually, it was hoped that instanton physics [^] might pave the way for a complete 
understanding here. For instance if we knew that the real gauge potential in nature 
was just the one-instanton gauge potential, and we knew how to calculate the density 
of states of the instantons, then we could calculate this-or-that condensate, on top 
of having direct information about the background field of vacuum fluctuations. By 
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studying the lifetime of the D meson, we will, at least, learn a few things about the 
character of these fluctuations. The D meson turns out to be an interesting system 
for the following reason: the charm quark is both heavy and light, so to speak. Heavy 
because a s (m c ) << 1 and light since in the OPE expansion, jj,/m c is not small, i.e. m c 
is heavy enough for us to use the OPE expansion, \ijm c ~ 0.5, but it is light enough 
that this expansion is not so good! Indeed, at the end of the day, we will find that with 
the constraints imposed on us for m c , and other pertinent parameters, the standard 
OPE treatment of the D lifetime FAILS! This failure of duality will then be analyzed. 
In the end, we will find that the D meson is not such a good kitchen for making the 
dinner of CP Violation, but is a good kitchen for making the breakfast of QCD. So 
what I will do in this thesis is little more than the standard QCD treatment - with a 
few bells and whistles which I helped develop over the past few years - of a heavy quark 
decaying with a light spectator. In the introduction, I will discuss the OPE method, 
reviewing how to calculate with it, and explaining why it is so useful for QCD, and in 
particular the decay of a heavy meson. Then, in the following sections, I will discuss 
perturbative and nonperturbative aspects of the problem of the lifetime of the D meson, 
and all revelant features and numerics thereof. In the end, we will be faced with the 
following conclusion: in the semileptonic decay of the D meson, duality fails, i.e. theory 
and experiment simply do not match. This failure should be taken as a warning sign 
for all other hadronic processes below ~ 2 GeV. For example, we can not trust the low 
energy a s extractions from r decays. 
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Introduction - the OPE 

"So, Dave ... you decided to study theoretical particle physics? Its very, very difficult, 

you know...." - Arkady Vainshtein. 

In the treatment of the decay of the D meson, we rely almost exclusively on Wilson's 
OPE ||. Wilson's OPE is actually nothing more than a method of calculating with 
a quantum field theory. It is just the method by which a non-local product of fields 
is expanded into a product of local fields. The utility of the OPE in QCD is almost 
immediately evident. The OPE provides us with a basis for dividing processes with soft 
and hard momenta into noncalculable and calculable parts. In QCD it is just what we 
need! The running of the strong coupling constant allows us to use perturbation theory 
for ultraviolet parts of diagrams, while we hide our ignorance of infrared dynamics in 
condensates. 

To learn the essentials of the OPE, it is sufficient to consider a simple problem in a 
simple theory. Everyone and his brother knows how to calculate the Green's function 
in scalar c/> 4 theory. Below, we will calculate it, and its first correction, with normal 
perturbation theory, side by side with the same procedure using Wilson's OPE. The 
example is stolen straight out of ref. fTpf . Basically, Wilson's OPE is just the following 
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D(q) 



x- 




Figure 2.1: The propagator for c/> 4 theory to O(A). 



relation: 



. oo 

/ e^dxT[j A (x)j B (0)} -+ ]T D£ B (q)O n (0). (2.1) 

n=l 



In the example below, we will implement the relation, see how it generalizes the normal 
perturbation methods, and directly see its beautiful utility in applications to strong 
coupled theories. For these illustrative purposes, we will consider a simple example: 
the Higgs model. We can assume the mass squared parameter of the Lagrangian has 
positive sign - there are no nonperturbative dynamics - and write the Lagrangian as: 



£ = jWf-^-^ (2-2) 



where mo and Ao are the bare mass and coupling constants, and the field (j> is real. Our 
goal is to calculate the Greens function, D(q), of the field <p in the one- loop approxi- 
mation using traditional methods, D pert and the OPE, D OPE . 

Consider figure 2.1. Using the Pauli-Villars regularization method, D(q) can be 
written as, 

D ^ iq) = -l^ [1 + ^ 2 -<^f/<\ (2 . 3 ) 
q z — thq 32tt z q 2 — rag 

Our task is to reproduce this result using the OPE. Essentially, the procedure is the 

same as the traditional calculation, except that the perturbative integrals - now the 



coefficients of given operators - are integrated down to some infrared scale jx, i.e. the 
perturbative graphs containing lines with momenta less then [i flowing through them 
are 'cut'. These cut graphs are nothing else than the operators. This particular example 
is an instructive one, since we can calculate the operators explicitly (unlike in QCD). 
The expression from the OPE we will need for correspondance is: 

D OPE (q) =< 0|Ci(/i)l + C^ 2 (/i)|0 > (2.4) 

Here, and in normal applications, the operator C\ stands for all of perturbation theory 
- with the caveat that we exclude infrared contributions. The excluded infrared contri- 
bution for C\ rears its head in the subsequent operator terms of the OPE expansion, 
each of which has its own perturbative contribution, ad infinitum. Note another key 
feature of the OPE - one which we utilize in heavy quark applications. The series of 
operators is one of increasing dimension, thus our coefficients must come along with 
decreasing dimension. The utility in QCD applications is that the expansion will be 
in powers ~ ^-qcd/Q where Q is some large external momenta. To wit, we expand 
our previous result, eq. |2.3| , for large q and thus have a prepared result in normal 
perturbative methods we can compare to our upcoming OPE calculation. Expanding 

for q 2 = —Q 2 , Q 2 >> m 2 , we get: 

11 \ M 2 

D per \q) - 7^2 + 7u \< + ^( m2 ~ < ln — )]■ ( 2 - 5 ) 



Q 2 Q 4 u 32tt 2 



rn, 



o 



Now lets get the OPE result. As stated above, the coefficient C\ is nothing but the 
perturbative result with the infrared cutoff fi - note here that the ease with which we 
implement the cutoff \i in general doesn't exist, for example a prescription for cutting 
off infrared dynamics for multi-loop processes would require great care (but is in fact 



possible [11 1). In the case here, the result is easily achieved - just subtract from the 



previous result, eq. |2.5| the same expression with the ultraviolet regulator M replaced 
by /i. Then, C\ is just 

11 A M 2 

Ci= -i? + w H + 5^ (M2 " " 2 - m ° V )] (2 ' 6) 



(a) 



(b) 





Figure 2.2: Figures corresponding to the OPE fi cutoff prescription for the calculation 
of the propagator in (J) 4 theory. 

Now we need to calculate the remaining piece. The coefficient, C^p, is just figure 2.2a. 
Calculating the diagram of figure 2.2a yields the result 



C, 



An 



(2.7) 



2Q 4 ' 

where the factor two comes from combinatorics. Then all we need is the vacuum ex- 
pectation value < (j? >. It is just the bubble graph of figure 2.2b with the entire 
contribution coming from perturbation theory. Note, here we remember to only inte- 
grate in the soft region from to fi, and so, 

<0|^ 2 |0>=^(/x 2 -mgln^|), (2.8) 

Combining these 'OPE' equations, we see that they give the same result as the standard 
perturbative methods, i.e. 



D 



OPE 



jypert 



(2.9) 



Turning to QCD, one quickly sees the utility of the OPE. The OPE allows us to use 
perturbation theory in the UV - where a s is small - and get some information from 
perturbative coefficients, and in addition, allows us to at least give our ignorance of 
the soft dynamics a name. Here we were able to calculate the operators directly. In a 
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strongly coupled theory like QCD, we have no idea what the Greens functions of quarks 
and gluons are for soft momenta and thus have no right to calculate them - instead 
we extract them phenomenologically, as was done first in SVZ sum rule applications 
]?]]. To be a bit more exact about actual practice, what is done in QCD is not exactly 
what is done in the example above. Technically, the coefficient functions - for infrared 
safe quantities - are integrated down to \x = 0. In the case of the D meson, it is just 
because we will borrow the the perturbative results of QED (muon decay). This is 
a blatent double count, but as it turns out a justified numerical approximation that 
works well in QCD applications. For instance, in the case of the gluon condensate, the 
true nonperturbative phenomenologically extracted result is numerically much greater 
than the perturbatively calculated contribution - a sign of the strongly fluctuating 



vacuum fields. This assumption goes under the name of the practical OPE [10]. The 
actual, proper, implementation of the OPE would result in the infrared contribution of 
coefficient functions being excluded, and thus the introduction of terms {[i/Q) n (see eq. 



2.6). The practical OPE, then, is the approximation wherein these terms are neglected 
- an approximation to be checked in each and every application. 



Chapter 3 

The Semileptonic D lifetime - 
r(D) - an Introduction 



"It is so easy to calculate it that it is impossible to make a mistake" - Misha Shifman 

If the D meson was just a point particle, all c quark, so to speak, then its decay 
wouldn't be all that interesting. Aside from its CKM matrix element, in fact, it would 



be the same calculation as for the decay of the muon, now a textbook favorite, e.g. [12]. 
Indeed, the semileptonic decay of the muon yields the decay width: 

IV - = (1 + A^a e ) (3.1) 

where the first order perturbative correction, A^ l \ was first calculated in ref. 13]. 

The calculation for the free quark decay c — > slv goes in the same way, since the 
weak currents governing the decay have the same form. Thus by multiplying the above 
expression by the CKM element squared, |V^| 2 , we have the D decay width with one 
loop corrections. We can play a quick numerical game and see whether the expression 
agrees with experiment. Setting a s = for now, 

Y(D) = & f\ V <»\ 2 ™* = L1 x 1Q -13 GeV (3 2) 

1927T d 
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(witb[]) m c = 1.4 GeV which is to be compared to the experimental value: 

^expt(D) = 1.06 x 1(T 13 GeV. (3.3) 

It looks like a great agreement! We shouldn't be so hasty, however. Corrections to the 
width will come both from perturbative and nonperturbative sources. Unfortunately, 
as we will see, the result gets spoiled by these corrections []. 

To calculate the nonperturbative corrections, it will be necessary to introduce an 
additional formalism, first proposed by Schwinger in QED [O], and later discovered 
independantly in QCD in |l(| with the first applications appearing in [17] and p8p . 



For now we focus on the perturbative corrections. the coefficient of the 0(a s ) 

corrections can be computed by considering diagrams of the type in fig 5.1 with gluon 
brehmsstrahlung occuring off of the heavy quark line or the intermediate quark. The 
result is 

* 1> ~SF<" s -f> (35 > 

If we wish to borrow this result for the QCD decay, we must address two issues. First, 
since the result was first calculated in QED, there was no reason to introduce the scale 
\x separating hard and soft effects, thus we need some theoretical justification for using 
the QED result in QCD where there is an infrared problem. Second, numerically it is 
an important issue as to what point we normalize a s . 

Formally, using the QED result, we imply the seperation scale \x = 0. Thus using 
the result in QCD, it means literally that everything to 0(a s ) has been calculated - 



1 This value of m c is taken from the early sum rule estimates I will return later to a full discussion 
of the heavy quark mass. 

2 As an aside, note that there is a phase space suppression coming from the fact that the s quark 
mass in the final state is massive. The suppression to the partonic width is given by 

<t>[z] = 1 - 8z + 8z 3 - z 4 + 12z 2 ln[z] (3.4) 

Numerically, taking the s quark mass as around 150 MeV, this effect is roughly 5 percent, if we are 
after an accuracy greater than this, then, of course, these phase space suppressions must be taken into 
account. 
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our work is done. Of course, this is nonsense - we have no right to believe that the 
quark/gluon Green's functions at low virtuality correspond to the strongly coupled 
dynamics that actually occur for scales less than Aqqd- On the other hand, we don't 
want to throw away this useful result. The answer to the puzzle was discussed in section 
2: we use what is known as the practical OPE |nj. Proper use of the OPE requires that 
we cut our integrals off at the infrared scale /x. Doing so would introduce corrections 



of the order (fJ-/Q) n (see eq.(2.6)). In the practical OPE, these terms are discarded. 
The assumption is that the perturbative piece that we have no right to include (from 
virtualities to fx) is numerically much smaller than the phenomenologically extracted 
operator. 

As for the second question, what scale do we use for a s , there seem to be two options. 
Formally, we are allowed to take whatever scale of a s we like since the difference is in 
the next order of a s . By choosing the appropriate scale, however, we minimize the 
corrections at next order. The most physically reasonable scale is obviously the scale 
fi = const, x m c . This is the scale involved in the physical decay - any emitted gluon 
would carry this momenta. To see which scale actually comes into play one has to 
consider next order corrections. Here, one can use what is called the BLM approach 



[19]. This generalized approach is applicable when one deals with a single gluon line 
(dressed by all bubbles). It amounts to inserting the unexpanded expression for the 
running coupling constant 

inside the integrand of a one-loop Feynman graph which depends on the gluon momen- 
tum k, with the subsequent integration over k. Letting a s run we have perturbative 
integrals like 

~/2£). (3.7) 

For the kinematics of the problem at hand, such integrals will be saturated at scales 
m c . 
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Now we would like to include nonperturbative effects. The most efficient way to 
calculate the nonperturbative corrections is by using the background field technique 
and the Fock-Schwinger gauge. In the next section, I review both formalisms. 



Chapter 4 



Calculations in external fields in 
QCD 

"Getting the damn 2's and it's in the right place is the whole point!" - Feynman 

Here, we review a formalism proposed by Schwinger j[4|] for the description of motion 
of particles in external fields. I will rely heavily on the formalism (developed in j^0|) to 
calculate the nonperturbative corrections to the semileptonic D lifetime. The Schwinger 
operator approach is formulated in coordinate space. Let us introduce a set of states 
\x > which are the eigenvalues of the coordinate operator X^: 

X^x >= Xf,\x > . (4.1) 

We also introduce the momentum operator P^ which satisfies the commutation rela- 
tions: 

[P fi ,X v ]=ig fll/ (4.2) 
[P^Pu] =igT a G% (4.3) 

where T a is the color group generator, related to the Gell-Mann matrices by T a = X a /2. 
In addition to the above relations, we will work in a specific gauge - the Fock-Schwinger 

15 
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gauge. The Fock-Schwinger gauge is the condition 

x^ = 0. (4.4) 

There are two conveniences resulting from this gauge choice that we will take advantage 
of. The first is that the potential is written in terms of the field strength tensor, G^, 
and thus easily yields gauge invariant expressions. The second is that A^(0) = 0, as we 
will now see. 

4.1 The gauge potential in the Fock-Schwinger gauge 

With our choice of gauge in hand, lets write out an expansion for the potential, A^. 
We start with the identity: 



A°(x) = / adaG^(ax) Xl> . (4.5) 



We can prove it, since, from our gauge condition, 

then inserting the definition of the field strength tensor, 
dA (y) 

Vp ~dy — = VpG w + y P d P A » + Vp\ A ^ A p\ = Vp G tJ.p + Vp d P A pi ( 4 - 7 ) 



so then 



My)+Vp*¥^ = VpG„ l {y). (4.8) 
Now, reparametrize y = ax, then 

^(a^(as) = Afj,(ax) + a-^A^ax), (4.9) 

and so, 

/ ax p GpJax)da = / — (aA^(ax)da = x^A^. (4.10) 
jo Jo da 
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Now, to get the expansion for the gauge potential, write the field strength tensor as, 
G^(ax) = G pp (0) + ax^^G Pfl (0) + ... (4.11) 
insert this into the integral representation and note that since in this gauge, d a = D a 
A°(x) = -La^Gp^O) + ^x a x p (D a G pp (0)) + .... (4.12) 

4.2 Calculating propagators in the Fock-Schwinger gauge 

As an exercise in this gauge, we can try and calculate propagators using the external 



field method [17|. We will need such propagators later. The scalar propagator for a 



massless particle in the Fock Schwinger gauge is written as, 

S(q) = J dxe iqx < x\-^\0 >= J dx< x\ —^ |0 >, (4.13) 



since 

e * x P ll = {Pr + qJe i « B . (4.14) 

Now expand this propagator in powers of P/q. Here we use a simplification resulting 
in our choice of the Fock-Schwinger gauge referred to above: 

yl M |0>=0. (4.15) 

Thus, our main strategy is to expand in P/q and pull all gauge potentials to the right. 
Also, since 

J dx < x\p p ■ ■■ = J dxdyid p 5(x — y) < y\ ■ ■ ■ |0 >= (4.16) 

we pull all p's to the left. The only other ingredient is just commutators between p's 
and A's, which are trivial to evaluate in the Fock-Schwinger gauge. Expanding equation 
O in P/q, 



1 _ 1 2Pq | 4(Pg) 2 P 2 | P 2 (2Pq) + (2Pq)P 2 8(Pg) 3 | Q(p4) ^ 



(P + q) 2 q 2 q A q & q 4 q e 
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and performing the procedure described above, we get 



D(q) = ^-^D a G ap q p 



~~z—g{q.aG ap Gpryq~ 1 + —q G ap G ap ) ^q~ i D~ i D a G ap q p + ... (4.18) 

Zq 4 q 



The case of the scalar propagator is of course easiest since we do not have to worry 
about the anticommuting 7 matrices. We will, however, need the result for spinor 
particles. The result^ (in the massless case) of this calculation is, to O(DG) : 



q 2q* 
3q 6 

3ig 7 jD 7 g a ,G f a p7 j£ ,7 5 ]. (4.19) 



-^-g [q 2 D a G ap j p - qD a G ap q p - q 7 D^q a G ap 'y p 



1 Note, here I use p instead of the Feynman convention, j/, to denote the Dirac matrices - Misha 
Voloshin taught me this, explaining that to him, a slash meant that the given variable was crossed out! 



Chapter 5 

Nonperturbative corrections to 
the lifetime of the D 

"The heavy quark is heavy - its like a locomotive flying through mosquitoes" - Misha 

Voloshin 

In section 3, we used the textbook calculation for muon decay to get our hands dirty 
for the case of the D meson. Let me derive this result in a fairly non-standard but 
equivalent way by using the optical theorem. The technique will be a bit long-winded 
for the standard partonic calculation, but will simplify life greatly when we are after 
nonperturbative corrections. 



5.1 Derivation of parton result using the optical theorem 

The optical theorem relates the imaginary part of a forward scattering amplitude to its 
observable cut process. The optical relation valid for the total widths is, 

r = -i-Im < D\f \D >, (5.1) 
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Figure 5.1: The forward scattering diagram for semileptonic decay. THe dashed line 
represents the 'cut', i.e. the imaginary part of the diagram 

(the relation above implies the use of relativistic normalization of states) which is our 
starting point. Here T is the time ordered product of, in this case, the two weak currents 
governing the semileptonic decay. From figure 5.1, we get the transition operator in 
configuration space as, 



T = ~ I d 4 x(5(r)i: /( .S' (/ (.r.())/ / ,,,(.r).[.VQ(0). 



(5.2) 



where Q represents the heavy quark external line, S q is the light quark propagator, the 
T's are the usual V x A Dirac matrices (r^ = 7 M (1 — 75)), and the lepton tensor is, 

2 1 



{2x»x u -x 2 g^). 



^27r 2 x 4 ^27r 2 a; 4J vr 4 x 8 
To get the parton result, take the light quark propagator as free: 

S q (x,0) = 



(5.3) 



(5.4) 



Upon convolution of Dirac matrices, T becomes: 

f = J d 4 xQ(x)f(x)(l+ l5 )Q(0), (5.5) 

where f{x) is a function of x. The interaction of the heavy quark with the light degrees 
of freedom enters through the Dirac operator, = — igA^t a . The background 
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gluon field is small compared to rriQ, and thus there is a large 'mechanical' part in 
the x-dependance of Q(x) p3|]: 

Q{x) = e im ^Q{x). (5.6) 

Our transition operator is still a nonlocal object. To write it as a local object we should 
Taylor expand the heavy quark field in the Fock-Schwinger gauge: 

dQ(x) = (1 + X a \ d a id a2 + ...)Q(0) 

= (1 + X a lD ai + ^X ai X a 2D al D a 2 + ...)Q(0). (5.7) 

The above expansion is legal since in the Fock-Schwinger gauge ordinary derivatives 
may be substituted by full derivatives at the origin. We can get rid of the 75 term in 
equation fl5.5|) since, through the equations of motion it appears as a full derivative 

< D\-^(Q W Q)\D >= 2im Q < D\Qj 5 Q\D >, (5.8) 

and the derivative of the forward scattering matrix element gives zero (this matrix 
element and others like it containing 75 are zero from parity arguments anyway). Upon 
convolution of the lepton tensor with T^'s, etc., we are left with 

f = const, j d A xe ipx Q(0)-^Q(0). (5.9) 

To do the integral (see the Appendix in pOfl), we use the result 

! =J^=Www iP ' ln[V1 ' <5 ' 10) 

which gives, in our case, 

~ d T . d i(-l) 5 2- 6 vrV(ln|p 2 | +m) , r „, 

r =-W = " S7 "* — m (5 ' u) 

and so finally^ we get 

L ~ l 384^ • (5 ' 12) 



1 Note, this, the near end result for the parton decay width, is the beginning of where we look for 
the nonperturbative 1/mj? pieces. 
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To leading order, p 4 = P 4 (we will see it later). Substituting this result, and using the 
Dirac equation, we are left with 

To{D) = A^ Im( ^S^ ) < D \Q(PM°)\ D > ( 5 - 13 ) 



Also to leading order (we will see it explicitly later (equation ( 5. 16j ) ) , 



< D\Q(0)Q(0)\D >~< D\Q(0) l0 Q(0)\D >= 2M D , (5.14) 

i.e. this operator merely counts the number of c quarks in the D meson - it is basically 
the number density. Here we get into points concerning fi dependence of operators 
involved in our expression. Technically, QQ carries (i dependence. Moreover, if // is 
large, say a few GeV, then we are not entitled to say that there is only one c quark in 
the D - we must evolve the operator down to a few hundred MeV - a typical hadronic 
scale. Later we will see that the heavy quark mass should be at a normalization point 
[i = m c . If so, then the operator QQ is also at \x = m c . We must evolve this operator 
down to a lower normalization point where we can say that all fluctuations are soft, and 
thus do not involve the creation of heavy quark anti-heavy quark pairs. Evolving this 
operator down is equivalent to pumping the fluctuation contained in the operator into 
the coefficient function. These fluctuations are actually already contained in the per- 
turbative correction, anyway, so we can forget about them. Now, upon combining 
the relations above, we finally get the expected result: 

5.2 0(l/mg) corrections 

To get the 0(l/mg) corrections, we will need one additional piece of machinery, the 
expansion 

QQ = Q^foQ — — ~~o~~2~~ ~ + total derivatives. (5.16) 
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2 



The 0(l/mn) result in the above equation was first derived in [21]. The relation itself 
can be easily established from the equations of motion. Since, 

^ C = i" C ' (5 - 17) 



then, 



which implies that 



-1-70 -1- 7o 1-70 1 r K1B v 
c^c = c— — c = ^cppc, (5.18) 



d{1 ~ 7 ° )C = + \° G)C - (5 - 19) 

Then, using the equations of motion again, we see that CTqC actually reduces to an 
operator of dimension 7, since 

and so we arrive at the result above. 

There are various possible sources for I/tuq corrections to the parton result. First, 
recall that the first term of the gauge potential in the Fock-Schwinger gauge is ~ G. 
This term actually gives a contribution of 0(l/mg) (we will see these terms later). 
There are however two other sources of 1/m 2 terms. The first is from equation ( |5.1£| ) 
above in the expansion of QQ, and the second is from the identity 

P 2 = p 2 - ~ a G (5.21) 



which we plug into eq. 5.12j (the identity above is obvious since 7^71/ = + o^ u ). 



Using these two results, we then plug T into the optical relation. The result is 

r 2 1/ 2 •? ,,2 1 ,,2 

T(D) = x [1 + A^a s - - ±% (5.22) 

v ' 192vr 3 2 m 2 2m 2 J y J 

Two points are in order here. First, notice that there is no contribution coming from 
operators 0(l/m c ), i.e. there is no operator of dimension 4 fl22|| . Second, since T is 
a Lorentz scalar, we should not be surprised to see that the only operators appearing 
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in r are Lorentz scalars (/^ enters only in combination with Q^foQ - from the original 
Lorentz invariant operator QQ). 

Our task is now somewhat complete - we have the partonic width, and the leading 
perturbative and nonperturbative corrections. Now we would like to do some numerics 
and see whether we match experiment. This is an area of hot debate for sure, but as 
we will see, the failure of the perturbative and nonperturbative corrections to give the 
experimental result is fairly insensitive to most of the numerical debates. Clearly the 
most important numerical issue will be the mass of the c quark, since it appears as the 
fifth power in the parton result. We now turn to a complete discussion of numerics. 



Chapter 6 



Theoretical meaning and 
numerics of nonperturbative 
parameters 

"...the set of values to be passed from the elders to the young generations included the 
idea that high energy physics is an experimental science that must be very closely 
related to phenomena taking place in nature...." - Misha Shifman 

"So Dave, you would like a problem to work on? Well, I must warn you that I will not 

give you something very mathematical, but instead more phenomenological - any 
results you will obtain will probably have direct relevance to current experiments..." - 

Misha Shifman, much to my excitement! 

In this section, I discuss the meaning and numerical values of the heavy quark mass 
m c , nl, and \j? g . 
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6.1 The meaning of the heavy quark mass 
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The issue of the heavy quark mass can easily be anticipated to be a sticky one. As 
everyone knows, quarks are not color singlets and thus do not appear isolated in nature 
(quarks are not asymptotic states, if you like). And yet, the theory of QCD is built 
upon quarks and gluons. Clearly, a key parameter in any QCD calculation will then be 
the quark mass - it enters in the QCD Lagrangian. In the calculation of the D meson 
lifetime, we will be borrowing perturbative results from QED: the one loop correction 
to muon decay. In the QED calculation, the muon mass used is just the pole mass. 
In QCD, this becomes a problem since by using the pole mass of the heavy quark, 
we are unjustifiably including long distance dynamics. Using the pole mass, we are 
essentially assuming that we have knowledge of gluonic Greens functions at all scales 

- a clearly unjustified assumption. In fact, it can be shown that the pole mass of the 
heavy quark is defined through an asymptotic series in a s , leading to an ambiguity of 
order Aqcd- This is a strong statement, and with the propagation of the words 'pole 
mass' in almost every paper on QCD, one gets the feeling that this fact is somehow 
avoided, or misunderstood. The misunderstanding comes from whether we are talking 
about the 'pole mass' or the 'one loop pole mass'. Although, as we will see below, the 
pole mass is ill-defined, one can rightfully use the one loop pole mass, as I will below, 
in any typical problem. In the problem of the D meson decay, use of the one-loop pole 
mass really helps simplify life - we can just use the well-known QED result for the muon 
decay, without including in it some scale dependance (i.e. cut off loop integrations at 
some scale). Then, use of the pole mass is proper, since we have, in effect, set the scale 

- at zero virtuality - by using the QED result. Here, we see the conceptual difference 
between the <p A exercise offered earlier, and what is called the practical OPE by QCD 
practitioners. Let me then, in turn, review two main points concerning the heavy quark 
mass: 
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Figure 6.1: The diagram corresponding to the mass renormalization due to the renor- 
malon gluon bubble chain. 

6.1.1 Illegetimacy of the pole mass 

Although commonly used in practical applications, it should be understood that the 
pole mass of the quark is an ill-defined object. The subject is thoroughly discussed in 



[26], and I only briefly touch upon some essentials here. To see that the pole-mass is 



ill-defined, consider the relationship between the pole mass and the running mass, 

™St = ™eWECn(^)(^r (6.i) 

71=0 

The above series is not well behaved, diverging factorially. Such a series is an asymptotic 
series, and must be truncated at a critical order in the expansion. After truncation, 
the remainder of the series shows up as an uncertainty (knowing the exact remainder 
would constitute having knowledge of the vacuum fluctuations). 

One particular source of the factorial divergence is the infrared renormalon. The 
infrared renormalon is the effect of the insertion of the bubble chain in gluon lines, i.e. 
allowing the running of the strong coupling constant, a s — ► a s (k 2 ), see figure 6.1. For 
nonrelativistic momenta, k 2 « fj, 2 , the expression for the mass correction is 

f d^k 4ira s (-k 2 ) f d 3 k a s (k 2 ) 
6m Q ~ / 7^73171 U / 732^7-- ( 6 - 2 ) 



{2-K)Hk k 2 J 4vr 2 p 

Expressing the running a s (k 2 ) in terms of a s (// 2 ), k 2 < fi 2 , 

a s (k 2 ) = aV)(l + ln^)- 1 ,b=^N c - \n } (6.3) 

4-7T n z 3 3 
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one can expand a s (k 2 ) in a power series of a s ([i 2 ). One then finds for the (n+l)-th 
order contribution, 

5m Q +1) i 2, u bajj^), n .... 

— — a s (n )7rn!( — — — ) . (6.4) 

m,Q An 

Observe that the coefficients grow factorially and contribute with the same signF]. There- 
fore one cannot define the sum unambiguously. An optimal truncation leaves one with 
an irreducible error of O(Aqcd)- To see this, note that truncating the factorial growing 
series at optimal order (see the section on duality violations in chapter 10) leaves an 
exponential uncertainty, 

A ("iQ°' e - ™qOm>)) ~ Mo ex P(~ faT7^y ) ~ A QCD- (6.5) 
6.1.2 Irrelevance of the pole mass 

Anyway, in applications, it is the running mass and not the pole mass which enters 
into our theoretical expressions. As explained above, the perturbative series for m P Q le 
diverges factorially. In our expression for the width, a similar divergence occurs in 
the a s expansion due to renormalons. If the running mass is used, however, both 
effects combine to cancel each other. This point is investigated thoroughly in |2q |. The 
result is really neat. It tells us, basically, to perform the following procedure when 



extracting quantities like the mass from total widths [33]. First, write the total width 
as an expansion in a s and l/m c . The a s series is, as in this case, to be taken at 
fj, = 0, i.e. the QED results can be borrowed. Next, one can numerically investigate 
the coefficients of the a s series. In the BLM procedure, the corrections at two-loops 
are easily calculated, and give large coefficients. Ref. tells us that these coefficients 
are nothing more than the signal of the factorially diverging perturbation series. Not 
to worry, since, by inserting the running mass, mq{mQ) we can 'cancel' the infrared 
renormalon effect responsible for these large coefficients. After this procedure, aside 



1 Same sign series can not be summed since the pole of such a series will lie on the real axis, and the 
contour around such a pole when Borel summing cannot be unambiguously defined. 
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Figure 6.2: The diagram corresponding to the running mass of the heavy quark. 



Figure 6.3: The diagram corresponding to e + e annhilation. 

from issues of duality violation, it is perfectly alright to extract the running mass from 
the theoretical expression for the width. The running mass itself is easily calculated 
(see figure 6.2): to one loop the result is 

2a 

™>q((j,) = m Q (no) + -ip^ ~f^- ( 6 - 6 ) 
6.2 The value of m c 

In this subsection we address the question 'what is the value of m c \ m c was first 
extracted from the QCD sum rule application to charmonium, and so we review the 
sum rules here. As a cross check, I also present an estimate of m c from mj. 

6.2.1 m c from sum rules 

Here, I derive one of the first OPE type applications in QCD - charmonium sum rules, 
|7]]. The derivation relies on the standard OPE treatment found in sum rules, or heavy 
quark applications. The starting point is quark-hadron duality: on the left-hand side of 
the equation, there is a theoretical construction in terms of quarks and gluons, on the 
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right hand side is a quantity extracted from some hadronic experiment. The assumption 
of quark-hadron duality is built in from the first step. As we will see later, this is a 
well justified assumption in applications where we rely on dispersion relations like the 
sum rules^j, but we must take care in making the assumption for the D meson. The 
reason for this is thoroughly discussed in chapter 10, and resides in the fact that sum 
rules rely on Euclidean kinematics, while the lifetime of the D meson is of Minkowskian 
kinematics. 

Our starting point in the investigation is the T product of two electromagnetic 
currents, (see figure 6.3): 

i < 0\J dxe^r^(x)i„(0)]|0 >= ( qft q v - q 2 g^MQ 2 ), (6.7) 

where = cry^c, and Q 2 = —q 2 . The polarization operator, II satisfies dispersion 
relations: 

d 1 f R c (s)ds 3sa c 

11 = / 71 i ^2 ' R c = T— 2" ( 6 - 8 ) 



dQ 2 12tt 2 Q 2 J (s + Q 2 ) 2 ' c 4vra 2 
We will especially be interested in the moments of II, 

M 1 f Rc(s)ds 1 d 

The simplest quark loop, is easily shown to be: 

(0) = J_ 2"(n+l)(n-l)! 1 

4vr 2 (2n + 3)!! (Am 2 ) n { ' 

Of course, the point of the sum rules was to go beyond such a simple expression, and 
include not only a s corrections, but also new power corrections like those we have been 
discussing above. In the case of charmonium, the first power correction is the gluon 
condensate. As it turns out, however, the contribution of the gluon condensate is highly 
suppressed for low moments, as are perturbative corrections - and so the key parameter 



It is another reason I review charmonium sum rules. 
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for low moments is just m c . Extracting it from the data, and ascribing a rather large 
uncertainty pi], 

m c (m c ) = 1.26 ± 0.10 GeV (6.11) 

6.2.2 m c from mj 

As a check of this result, we can examine another source where m c can be extracted. 
In |24|] the following result was first derived, 

M B -M D = m b -m c + ^-lA _ tifA + (m 2 c ,m 2 b ). (6.12) 

The meson masses are measured accurately, (Mb = 5.287 GeV and Mjj = 1.864 GeV), 
and the value of mj, is known with high accuracy from upsilon sum rules [3^] mb(lGeV) = 
4.64 ± 0.05 GeV. Plugging in these values, we get m c = 1.2 ± 0.05 GeV where the un- 
certainty comes mainly from the uncertainty in [i 2 . (see below). 

6.2.3 Comments on the literature 

Despite the above arguments, there are still claims in the literature that one can take 
the pole mass of the charmed quark to be very high, m c = 1.65 GeV. A typical 
argument can go as follows. First, expressions for T(D), and mP° le are given, both with 
full one-loop and BLM type two-loop corrections. The explicit expression for the mass 
is 

m P^ = rn c [l + l^£l(l + 1.04)] (6.13) 

3 IT 

where '1.04' is due to the BLM type two-loop correction. The high value m c = 1.65 GeV 
is obtained by directly accomodating the experimental measurement of T(D) with the 
full theoretical expression, including nonperturbative effects. Then one uses the above 
equation to run mP° le to m c , and obtain m c ~ 1.34 GeV which is claimed to be close 
to the sum rule prediction. 
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A clear signal that this analysis is in trouble is that the value of the pole mass of 
the b quark extracted is m^ ole ~ 5.0 GeV in disagreement with Voloshin's extraction, 
m poie ^ ^ g Q e y What is wrong with the above analysis? The error is nothing more 
than a manifestation of the problems occuring when using the pole mass. Clearly, the 
above expansion for m c is not converging well. This fact is actually to the benefit of 
those who would like to fit mP° le from T(D) exp t and then run down to m c - the factor 
'1.04' allows such a large running. The problem is that it is not the pole-mass which 
really enters theoretical expressions, but the running mass. 

Above, it was discussed that the proper mass to use in expressions for the full width 



is the running mass. In ref. [ 33 1 it was explicitly shown that using the proper running 



mass, we reduce large second order corrections to the width in heavy quark decay. 



Explicitly, ref. [33] showed that for the semileptonic width we have large two-loop 
corrections, 

M = 1+ , ( S )+ ^). (6 .14) 

Yq[B) IT TT 

where a,2 is the BLM type correction, and a\ = —2.41, and 02 = —19.7 Q but upon 
insertion of the proper OPE scale, say, /i ~ 400 MeV , a\ — > —0.5 and 02 — * —0.4, and 
thus upon insertion of the running mass, the second order corrections become small. 
Upon insertion of these corrections to the total with, we obtain 

T(D) = r (l - 0.07 - 0.008 - 0.27 - 0.09) = r (0.55), (6.15) 

where —0.07, —0.008, —0.27, and — 0.09 are from a s , a^, /j,q, and corrections respec- 
tively (I will discuss these numerics in full in the next section). 

Now, extracting the mass from this result, the running mass mind you, we get 

m c (400MeV) = 1.57GeV (6.16) 



3 the result for ai is presented in the V scheme, a value around ~ —30 would be obtained for the 
MS scheme used in some papers, allowing an even GREATER value of m c . 
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which is in disagreement with the sum rules estimate, a result which I discuss again in 
section 7. When using the proper OPE procedure, outlined above, the extraction of m c 
from semileptonic D decay and charmonium sum rules disagree. 

6.3 About n 2 G 

The mesonic matrix elements of the chromomagnetic operator, oG can be extracted 
from the hyperfine splitting. First, note that 

B-G = B- B + 0{l/m c ). (6.17) 

Now, B is an axial vector and is determined by the dynamics of the light spectator 
cloud, thus we have only one choice, B = const x S q and so the matrix element must 
be: 

< Qvl^^-IQV >=C< Qq\S Q ■ Sg\Qq > . (6.18) 

Now, 

S~Q ■ Sg = \(Sq + S q f ~ \{S% + Sf ) = 1 -S TO t{Stot + 1) - \, (6.19) 

and, 

< (T|<7 • B\0~ >= -3 <l~\a- B\\- >, (6.20) 

since for an uncorrellated spin state, the matrix element is zero. Then, we find simply 
that, 

4 

M (0 -) - M (1 -) =< (T|<r • £|(T > - < \~\ff ■ B\\~ >= - < 0~\a ■ B\0~ > (6.21) 
and so, 

< i? G > D =< D\(Q l -a ■ G)Q\D >= ^m c (M D * - M D ) = -^{M 2 D , - M 2 D ). (6.22) 
Numerically, taking this mass splitting from the D and B system respectively, we get: 
< fa >D = 0.41 GeV 2 , < f x 2 j> B = 0.37 GeV 2 . (6.23) 
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A measure of the reliability of the nonperturbative expansion is then provided by 
\J < fj,Q >d /n^c = 0-46. And thus, right off the bat there is good reason to pursue 
higher order nonperturbative terms, 0(l/m^), in the expansion. 



6.4 About fj% 

The mesonic matrix element of the time dilation operator, < D\Q(iD) 2 Q\D >= is 
not known accurately. An analysis based on sum rules yields 



»l = 0.5 ±0.1 GeV 2 . (6.24) 



This result is in agreement by a bound from [27], [ 28 [ , |29], [34] 



/4 > /4 > 0.37 GeV 2 (6.25) 



This inequality is of huge importance, and much controversy [31]. The derivation of this 
inequality using quantum mechanics, taught to me by Misha Voloshin in his excellent 
class on phenomenology, is simple. Consider the absolute positive quantity, 

(a ■ vf) 2 = vf 2 - go ■ B, (6.26) 

then, obviously, the above inequality holds. The field-theoretic inequality was later 
derived in ||34| . There, the sum rule for the pseudoscalar weak current, J5 was written 
at zero-recoil, 

hlr^'^^^-^ <6 - 27) 

Since the structure function is non-negative, one arrives at the conclusion that 
/4(/^) > Mg(a*)- 

Still, there are arguments in the literature that this inequality is broken by the 
perturbative evolution of /U 2 . Consider figure 6.4 below The figure represents the issue 
of concern regarding the impact of the perturbative evolution on the numerical value 
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Figure 6.4: The diagram corresponding to the perturbative evolution of /U 2 . 

of fi"^. I bring up the issue not only due to its possible numerical relevance, but also 
because it sheds light on the actual meaning of both and rriQ, and their relation to 
each other. Recall that in the case of the operator QQ (i.e. just the heavy quark mass 
diagram) we had a result (for the running of itiq) which had no logarithmic dependance. 
In such a case, it is tempting to - due to the safe infrared limit - run fi to zero. For 
instance, this is what we have done in the case of QQ. As we have discussed above, 
however, this is an illegal procedure. As stated above, however, what is done in practice 
is the following: one takes known expressions, in our case one- loop QED corrections, 
and adds to these perturbative terms nonperturbative corrections expressed via certain 
matrix elements. What is done in the case at hand is the following. 

From figure 6.4 it is easy to see that we can define the one-loop perturbative ^ as 
just the infrared piece of the operator QQ, i.e. 

Q^ 2 Qone-loop = ^-I^QQ (6-28) 

Now, as we know, this quantity should be subtracted from the actual /x 2 (we cannot 
double count this piece!) Numerically, one usually considers a value of /U relatively high 
so that a s is small, but yet small enough so that the expansion in fi/m c is good. To 
estimate the perturbative evolution, I will choose /j, = 1 GeV. Doing so, the correction 
to is quite sizable - 0.15 GeV 2 . The inequality stated above, > y? G still is safe, 
however, as can be seen from eq. 6.27, at any scale. 
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For our purposes concerning the D lifetime, we realize, however, that the theoretical 
value of may not be completely crucial, for even if we take // 2 to have a very large 
perturbative evolution, so that /x 2 = 0.15 GeV 2 , and so assuming the inequality held, 
/Jq = 0.15 GeV 2 still we would have trouble fitting the lifetime, as can be seen in the 
next chapter. 



Chapter 7 

Putting it all together- The 
numerical prediction for T(D) 
with leading order perturbative 
and nonperturbative corrections 



With the issues of numerics under our belts, we can finally check to see how our formula 
for T D , 

r ( fl) = r„(i + A ia .-^-^), (7.D 

works. We take the mindset that we push everything in the direction of agreement. This 
corresponds to a high m c = 1.4 GeV, low n\ = 0.5 GeV 2 - 0.15 GeV 2 = 0.35 GeV 2 
and low /j,q = 0.35 GeV 2 - I slightly lower the value of ^ to reflect the fact that 
the inequality stated above still holds (due to the perturbative evolution mentioned in 
chapter 6). Plugging everything in to the lifetime, then, we get: 

T{D) = r (l - 0.24 - 0.27 - 0.09) = r (0.40) (7.2) 
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Recall that for m c = 1.4 GeV we had perfect agreement with the lifetime without 
corrections, thus the corrections spoil agreement, lowering the lifetime about by half. 
To alleviate such a situation, we have to really push our numbers: one choice is to 
make m c = 1.6 GeV (1.6/1.4) 5 « 2). This raises the rate by about half, and solves our 
problem. However, doing so really kills any agreement with the m c extraction from the 
QCD sum rules. Under our philosophy, we wanted to take the sum rule determination as 
bedrock - especially, since as we will see there are exponentially small duality violation 
possibilities for sum rules, but not for decay widths. 

Another option is to try and drastically lower ^ For instance, it might be that 
the extracted value of /x 2 = 0.35 GeV 2 then the perturbative evolution could knock the 
value down by say 0.15 GeV 2 , and assuming that the inequality still holds, we still get 
a large theoretical deficit: 

T(D) = r (1.00 - 0.24 - 0.153 - 0.05) = r (0.55) (7.3) 

And so it seems that even really pushing our nonperturbative corrections to their lowest 
possible values, we are still in trouble. 



Chapter 8 

\jvr? c Corrections 

"How tragic is wisdom when it brings no profit to the wise" 

Now, it looks like we are stuck. How can we accomodate such a large width? A 
fairly non-imaginative option is that for some reason the - it is always an option - 
perturbative, or nonperturbative corrections are enhanced at the next order. Thus, for 
instance, someone could try to go farther and calculate the \jrr? c corrections. This 
attempt was made in |3q| , and in this section we review the calculation of these cubic 
terms. As we saw in the last section in the numerical discussion, the hope is that the 
1/m^ corrections can give us an increase in the width around 30% to 50%. Clearly, 
such a contribution will be due to an unsuspected enhancement in the coefficient of the 
operator. 

Since the total width is a Lorentz scalar, the only new operators relevant at the level 
of dimension 6 are the four-fermion operators of the type: 

6 = cTqqTc (8.1) 
where q stands for the light quark, and T stands for Lorentz, and color structures. There 
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Figure 8.1: The four-fermion term in the transition operator as it appears at the level 



a 



are two distinct sources of the corrections, just as there were for the 1/m 2 correc- 

tions: operators of dimension 6 arising from the expansion of T, and l/m c corrections 
in the D meson matrix elements of the operators caGc and cc. 

8.1 The four-fermion operators at 0(a s ) and in LLA 

A four-fermion operator appears in T in the zeroeth order in a s in figure 8.1. The 



corresponding result can be read from Eq. (17c) of ref [37|, 



G 2 m 2 

Im f = [c^cfc - {2/3)c iW c k ][s k r^ i }, (8.2) 

here = 7^(1 — 75). The expression above, however, yields zero for two reasons. First, 
in the factorization approximation, T° corresponds to the annhilation contribution cs — > 
lu, and our spectator here is not an s quark. Second, even if we deal with the D s , the 
chiral structure yields a vanishing contribution when factorization is invoked, since 
under factorization < D\T°\D >= 0. Thus this is not the contribution that we are 
looking for. But lets not be hasty. The effective Lagrangian given above is at the 
virtuality scale, \x = m c . The type of pre-asymptotic effect we are after is at a different 
scale, fjL = Aqcd- Under the rules of the OPE, we therefore must 'evolve' the operator 
down to this infrared scale ('evolving' corresponds to sloshing some of the operator at 
scale m c into part of the coefficient at the scale Aqcd)- Physically, this corresponds to 
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Figure 8.2: The penguin graph for the four-fermion operator in T. 

including gluonic degrees of freedom, and thus, we can expect that the chirality, and 
hence vanishing character of the contribution which we obtained at scale m c will be 



different at the scale Kqcd- This contribution is calculated in ref. 37], and is 



Imf^-^lK/t-^ln^l^T^c+icf^c)^^^], (8.3) 

here f M = 7 M (1 +75). 

Notice that now, the light-quark current runs over all flavors and both left-handed 
AND right-handed fields appear, leading to the non-vanishing contribution of this result 
after factorization. This is a typical situation with the penguin contribution, and arises 
from the simple fact that the gluonic correction connects the V x A strange quark 
current to the pure vector background current. Above we have calculated the 
contribution arising from the logarithmic mixing of penguin operators and operators 
built from the light quark currents. Keeping only logarithms, however, may be an 
unjustified assumption, since m c //j, is not large. Thus, in the next sections, we calculate 
the full 0(a s ) contribution. 



8.2 Full 0(a s ) calculation 

Here we divide the calculation of the full 0(a s ) contribution into three seperate pieces. 
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8.2.1 The 'free' piece 

To start with, we calculate l/fUg pieces coming without any logarithms. We get terms 
here in the same we got 1/m 2 , terms when we examined Qp 4 pQ. And so, we can just 
start with the object: 

r m f - G 2 f \Vcs\ 2 Q(0)p 4 pQ(0) . . 

ImT °- 384^3 > ( 8 - 4 ) 

where 

= iD^ - gAp. (8.5) 

There are two possible contributions coming from the first two terms of the expansion 
for Ap, recall that 



(higher order terms will not contribute dimension 6 operators). First, lets focus on the 
firstQterm of the expansion of A^, and so consider 



QPP A Q = Q{P - A){P - AfQ 

= Q(p-^ Xp G m )(P-^x p G p ^Q. (8.7) 

As before, the strategy is to drag all A^s to the left, and P^s to the right. When Ap 
stands to the left it gives zero (Q is taken at the origin, so since A ~ x, AQ — > 0), and 
when P stands to the right, the Dirac equation can be invoked. Thus, we can already 
write 



Q(P)(P- -x p G ptl ) 4 Q. 



Continuing the procedure of pulling all gauge potentials to the left, we will find ourselves 
in need of a few identities: 

1 i 

[Pp, Aft] = iDpAp = iD^-XpGpp) = -XpDpGpp, (8.9) 



1 In re f. frS5fl this contribution was left out, but has minimal numerical effect anyway. The discrepancy 
between P5|]and |36| still exists even after this left out contribution is included 
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[P„, A„] = iA,^ = % -{G^ + x p (D v G plt )), (8.10) 

7/,G^P, = -^(PtG - <tGP). (8.11) 

Using these identities, the final result for the first term of the gauge potential expansion 
is 

Ql P D a G ap Q. (8.12) 

The second term in the expansion of the gauge potential also will give us contribu- 
tions 0(l/mg). Again, consider 

Qp 4 pQ = Q{P-Af{P-A)Q 

= Q(P 4 - P 2 P ■ A - P 2 A P-A- PP 2 - P ■ AP 2 )(P - A)Q (8.13) 

We drop the term ~ P 4 P since it gives no terms 0(l/mg), and also drop the term 
~ A ■ PP 2 since A here is completely to the left. Thus, there are three terms which we 
need to examine which will give 0(l/mg) contributions: 



Term 1 : Q(-P A A)Q 

Term 2 : Q{-P 2 [P^ A p ]P)Q 

Term 3 : Q(-2P 2 ,4 • PP)Q 

(8.14) 

To get the results of each term, one merely needs to know the few simple commutators, 
easily worked out above in the Fock-Schwinger gauge. The results are, term by term, 

Term 1 : ^Q^G^^Qm^ 
Term 2 : ~Q{D 7 G 7 ^)Qm 2 Q 

Term 3 : ^Q(D 7 G 7/i7/1 )Qm|. (8.15) 
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Adding the results from both the first and second terms of the expansion of the 
gauge potential we get the transition operator: 

f\ free = -iG 2 F \V cs \ 2 ^m 2 c c lp t a cq lp t a q, (8.16) 

where we have used the equation of motion, 

D a G a a ^ = -gq^t a q. (8.17) 



This piece contributes, upon factorization (see later in this section) 

Ar _ Aa s nflM D 



(8.18) 



T 3m 3 c 

Now we calculate the full logarithmic contribution referred to above. It will carry 
with it the piece (infrared divergent) proportional to log that we treat with the OPE 
cutoff prescription, and will also carry a finite piece. 

8.2.2 The 'Log' piece 

There is another contribution in the expansion of T arising from O(DG) terms in the 
expansion of the light quark propagator. The origin of the term proportional to D a G ap 
was made evident in section 4. Here we need to keep the mass of the s quark finite, 
however, since we will use it as an infrared regulator (we will perform the proper OPE, 
subtracting out from our coefficient the soft contribution form to fi). Keeping the 
mass of the s-quark finite results in the introduction of McDonald functions into the 
propagator. The result is: 



S q (x,0) = 



im 2 Ki(m\/—x 2 ) 



V—x 2 



m 2 x j , Gp A mK^mV-x 2 ) 

+ j^2mK (mV -x 2 )a pX 



+ -g^(2K (mV^x 2 )D a G ap7p - {D a G ap xx? 

+ xYfl 7 G a /-3 K Vi)to)^^)l- (8.19) 
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DG DG 

Figure 8.3: The diagram with the DG term in the s-quark line. 

Due to the introduction of the Mcdonalds functions, we will have integrals of the type 



e ipx d i 



X 



(8.20) 



(-x 2 )p 

which can be found in ref. Ipq] . 

To get the 0(l/m^) pieces, we isolate the piece of the propagator proportional to 
DG and one gamma matrice, see figure 8.2. The singular nature of the contribution 
can be understood, since K n (m s (V — x 2 )) ~ ln[m s (\/— x 2 )]. 

To avoid the infrared singularity, which shouldn't be included in the proper OPE 
calculation anyway, we calculate the graph of figure 8.3b, and subtract it from the graph 
of figure 8.3a. In this graph (10b), the lepton pair is hard - pointlike - and the s quark 
is soft. After subtracting the OPE infrared piece, we get the result, 



a 



= ^|^|(^m^(ln-f + -)) 
+ {2cT„t a c + cf^c) qi»t a q 



121) 



which gives a contribution aftor factorization, 



Ar 



9m 3 . 



(In^f + {)f D M D . 



(8.22) 



Lastly, we cannot forget to include pieces proportional to DG coming from the 
expansion of 1/m 2 operators. (Note, it is exactly how we picked up a piece of the 1/m 2 
contribution but there are no 1/m 3 terms coming from cc). 



46 

8.2.3 The '1/m 2 ' piece 

Just as there were contributions of 1/m 2 pieces coming from leading order terms so 
also are there 1/m 3 pieces originating from next to leading order terms. Previously, we 
expressed < D\c^oGc\D > in terms of the D*D mass splitting. This result is only valid 
to the leading order in l/m c . Let us observe that the spin splitting yielding M D » — M D 
is determined by the following terms in the heavy quark Hamiltonian: 

AH = —o ■ B + — ^r<? • E x 7? , (8.23) 

where B and E are the chromomagnetic and chromoelectric fields respectively (B = 
gB a t a and E = gE a t a .) To leading order 

A D = ^{Ml* - M D ) = - <a-B>= 0.405 GeV 2 . (8.24) 

At the level of l/m c the second term in the heavy quark Hamiltonian becomes important 
in Ad, as well as the second-order iteration in (2m c )~ 1 < a ■ B >. Assuming that both 
effects are of the same order of magnitude, we can roughly estimate the matrix element 
(2m" 1 ) < a ■ E x 7? > as the difference between A/j and A#: 

|(2m c )- 1 <<T-^X7r>|<A D -A i3 ~ 0.04 GeV. (8.25) 

Next observe that 



1 _ 1 1 
-ca uu Gn U c = -co • Be co • E x ire ciDiEAc. (8.26) 

2 m c 2m, 



The last term in the above equation reduces to the four-fermion operator which we 
take into account explicitly. The second term will be estimated as an uncertainty in 
the expression relating < ^oG > to Ap: 

=< \o G >= A D ± 2(A B - A D ) - (2m c )- 1 47ra s < c llx t a cq^t a q > . (8.27) 

Using factorization for the Oq term above, and the same values for the parameters as 
above, we get +0.01 for the contribution for Oq, so that 

H 2 G = A D ± 2(A B - A D ) = 0.42 ± 0.08 GeV 2 . (8.28) 
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As for /xj, it was shown in 2B] that the sign of the l/m c correction is negative. We 



will assume that the error bars in the numerical value of itself give the estimate of 
its l/m c contribution. 

8.3 Tallying the result 

In estimating the numerical value of the contribution we choose fr> = 170 GeV, 

a s = .3, and m c = 1.4GeV. Adding up all contributions, we get 

Ar 

— = -0.06 ±0.06 + 0.03 (8.29) 

here the first number is due to the four-quark terms in the transition operator, the 
second is due to the uncertainty of Oq, and the third is due to O n . Unfortunately, the 
result doesnt help our situation. Here we hardly have the 50 percent contribution we 
were looking for. 

So it seems that again we are stuck - is there no way out? 



Chapter 9 

Ways out 

"Deny Everything" - Hunter S. Thompson 

It appears that, in the case of the semileptonic decay of the D meson, all of our methods 
of attack have failed. What could have gone wrong? In this section, we review the case 
up to now, focusing on possible solutions. Some possible solutions were reviewed in 
previous sections, and relied on increasing or decreasing relevant numerical parameters. 
There, we found it is clear that we cannot increase or decrease the numerical value of 
any parameters in the expression for T(D) to obtain agreement with experiment - at 
least not at the cost of creating even deeper puzzles in other charm phenomena. Also, 
it seems clear that the way out of this problem is not due to the perturbative series. 
In previous sections we saw quite clearly that the next order perturbative corrections 
make minimal impact. What other possibilities are left? 

One approximation that was made in chapter 8 was the factorization approximation 
in the evaluation of the dimension 6 matrix elements. Corrections to the factorization 
approximation were reviewed in j39|. There, corrections were shown to be proportional 
to 1/N C and seem to be small. Remember, we would need for the corrections to be 
quite large. 

The next possibility might be that dimension 7 operators could save the day. In 
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principle, this may happen since the expansion parameter is ~ yjj^/m c ~ -5 and is of 
order unity. However, since the correction due to terms is of order 10 percent, 

this seems unlikely - not to mention, it would involve a very painstaking calculation! 

The most likely solution, however, does lie in the fact that our nonperturbative 
expansion parameter is large. This idea is tied conceptually to the fact that the kine- 
matics of the problem at hand are essentially Minkowskian. One justifies an OPE 
based procedure by keeping in mind an analytical continuation. In the problem of the 
semileptonic width this may be a continuation of the lepton pair - one considers the 
transition operator at such momenta where one is actually off the cuts corresponding 
to production of the hadronic states - in the Euclidean domain. The prediction on the 
cuts is made by invoking dispersion relations, in full analogy with what is usually done 
in the problem of the total e + e~ annhilation. In general, one can analytically continue 
in some auxiliary momenta which has nothing to do with any of the physical momenta. 

Whatever analytic continuation is done, the prediction for each given term in the 
l/m c expansion refers to the Euclidean domain and is translated to the Minkowski 
domain only in the sense of averaging which occurs automatically through the dispersion 
relations. If the integrand is smooth, however, we can forget about the averaging 
because in this case smearing is not needed. This is what happens, in particular, with 
the total hadronic cross section in e + e~ annhilation at high energies - quark-hadron 
duality sets in and the OPE-based consideration yields the value of the cross section 
at a given energy, locally (without smearing). At what energy release is the integrand 
smooth and can the terms in the l/m c expansion be predicted locally? The existing 
theory gives no answer to the question, but the problem at hand seems to suggest that 
the duality limit sets in well above 1.4 GeV 

In the next section, I review the work done in |4^] regarding duality violations. This 
work was the first real attempt to try and capture the essence of QCD violations, and 
provides probably the best framework with which to escape from the puzzle at hand. 



Chapter 10 



Duality violations 

"Dave - you can't just learn this subject by reading papers about it, you have to solve 
some problems by hand!" - Misha Shifman 

"You don't get it? Sit down with the book - AND A PENCIL IN HAND! - and work 
out every intermediate step, and you'll get it. Good for the soul you know...." - Serge 

Rudaz 

Depressing as our failure may seem at first, we should not look at this result as a 
failure, but instead as an opportunity. In fact, the D is an ideal testing ground for 
duality violations. Our QCD based calculations have failed to describe the hadronic 
dynamics, lending evidence to the idea that there must be a contribution which we have 
left out. In this section, I discuss the concept of duality, its violations, and its relation 
to heavy quark physics - in particular the problem of the semileptonic D decay. 

The concept of duality was introduced long ago j5| . Essentially, its meaning is this: 
we calculate amplitudes with quarks and gluons, and compare our results to hadronic 
observables. If duality holds, then the complicated hadronic dynamics involved in a 
given process yield essentially little effect - the most important effects come just from 
the simple point-like interaction of the quarks and gluons. In accounting, say, for 
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the effect of Fermi motion - including /x^, or in general any other operators, we take 
a stab at including some of these hadronization effects - from the quark-gluon side. 
Thus, each time we include some new operator, etc. into our calculation we redefine 
what we mean by duality. In the present problem of D decay, our definition of duality 
consists of the expansion in a s and l/m c corresponding to the hadronic lifetime. We 
prayed that this quark-gluon calculation was dual to the hadronic decay width. Well, 
it looks like its not! Thus, we need to stretch ourselves, and find some contribution 
which we have left out. Clearly it is a difficult task. Since we have the first few 
terms in both the perturbative and nonperturbative expansions, the effect that we are 
looking for will have to come from somewhere deeper. The perturbative expansion goes 
like ~ a s /ir which is numerically ~ 0.1. The nonperturbative series goes something 
like J < Hq > /ml ~ 0.5, so the nonperturbative series seems to be most suspect if 
our duality analysis targets poor convergence properties of our expansions. The first 
suggestion of this type of contribution came in ref. [f47|]. 

10.1 Exponential terms and asymptotic series 

Consider the series Yin a n n - This is an example of an asymptotic series: for a few first 
terms the series homes in on a limiting value, but then proceeds to skate away with 
the inclusion of higher order terms. Of course, including all the terms (for the problem 
at hand) we get a finite result. In the perturbative and nonperturbative expansions 
we deal with this type of series. Actually, it is not a scary thing that we deal with 
asymptotic series: hopefully, it can be the case in some situations that the impossible 
task of calculating higher order terms in a given expansion is indeed a groundless task. 
Perhaps in a given case the optimal order is n = 2 or 3, and we can forget about the 
higher order terms - the series should be truncated here. This might be just the case 
with the D meson - not a bad hunch, since we have experience leading us to believe 
l/m c is a large expansion parameter. 
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Truncating the series at a finite order, we introduce an exponential error. This can 
be seen in the following way. Suppose we have some function, f(x) represented as the 
expansion, 

f(x) = a + a\x + a 2 x 2 + .... (10-1) 

if the expansion is factorially growing, the coefficients can be written like a n ~ c n n\ (in 
QCD, c is roughly //). Now, as stated above, the series converges, but then starts to 
diverge again at an optimal n. Roughly this n occurs when 

x ~ (^n!) 1 /™. (10.2) 

We want to take into account the size of the last term neglected. Using Stirling's 
formula, 

n!~V2^n n e" n , (10.3) 

it is easy to see that the last term neglected is ~ exp(— x/c). Thus, truncating the 
series in fi/niQ we get an exponential accuracy ~ exp(— ptuq) where p ~ 1/Aqcd is 
some infrared distance. 

10.2 Physical picture of the exponential terms 

Surprisingly, there is a clear physical picture of the exponential error |fjjf| . Consider 
figure 10.1. In discussing the figures, let me reemphasize some of the points I have 
made earlier regarding the OPE. First, consider figure 10.1a. Figure 10.1a is a usual 
perturbative diagram where the quark emits and reabsorbs a gluon. As is typical of 
perturbative graphs, a large energy release is carried by one or at most two quanta. 
Figure 10.1b is the nonperturbative 'cut' diagram of figure 10.1a necessitated by the 
OPE. The cut line (cut since we have integrated down to our infrared point /i) now 
becomes an operator, in particular, p 2 . Figure 10.1b. represents the diagrams which 
we normally take into account, and have, in the case of the semileptonic D lifetime, 
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Figure 10.1: Possible diagrams contributing to heavy quark decay in a toy model. 

where one or at most two lines become soft. Figure 10.1c represents an additional 
contribution, the kind which gives us the exponential effect that we are looking for. 
Here, we have neither one or two hard lines nor one or two soft lines but rather, a 
hard contribution transmitted by many lines, so that each line is soft. The situation is 
one where a hard momenta is shared by a coherent, possibly classic field configuration. 
As we will see, by using instantons - just a solution to the classic equations of motion 
- we will be able to model the exponential contribution which is conceptually related 
to figure 10.1c. Actually, these type of exponential terms were known long ago since 
the early days of QCD f7|. There, they were essentially disregarded as most of the 
early QCD application, e.g. sum rules, are of Euclidean nature (dispersion relations 
are used). In typical heavy quark cases we deal with Minkowski type kinematics (no 
dispersion relations, but instead a direct analytic continuation to the physical cut), and 
here the nature of the truncation error differs drastically - the exponentially decaying 
truncation error in the Euclidean domain becomes oscillatory in the Minkowski domain, 
suppressed only by powers of our expansion parameter. 

Our task then is simply to develop some methods to generate this type of exponential 
term. In fact, the exponential terms appear in the same way as the power terms before 



them. A concrete example is given in [26]. There, it is shown that the truncation of 
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the factorially diverging series in a s leads to an exponential error 



Thus, from the exponential error involved in truncation of the perturbative series, 
we can see the presence of power terms. Here, we consider the series of power like 



terms. This too is an asymptotic series [41], and its truncation will give an exponential 



contribution of the type we are interested. 

10.3 Generating the exponential contribution 

To illustrate the main details of the instanton calculations, lets outline the practical 
motivation for the inclusion of the corresponding effects from the general perspective 
of the short distance expansion. Consider a generic two point function U(Q 2 ), the 
polarization operator of two vector currents (the Lorentz structure of the currents is 
irrelevant). 

II(Q) = J d 4 xe iqx Ii{x) = - J d 4 xe iqx < G{x,0)G(0,x) > (10.5) 

where G(x,y) is the quark Green's function in an external gauge field and averaging 
over the field configurations is implied - we work in Euclidean space. The normal power 
type corrections occur when we consider the expansion of Green's functions near x=0 
where the Green's function is singular. 

We now examine the question 'what happens when the propagator has a pole not 
near x=0?'. For example, consider the simplest finite-x singularity for the two-point 
function, 

U ( X ) = 7 27 2W ( 10 - 6 ) 



then, 



n(Q 2 ) = J d 4 xe iQx U(x) ~ e~ Qp . (10.7) 
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These types of contributions are not seen in the normal OPE, and thus represent a 
violation of duality. Moreover, they mimick the type of behavior we expect from the 
higher order terms left out in the truncated series. 

Unfortunately, the physics of these duality violating terms brings us, again, to the 
main difficulty of QCD - we would like to know the background field fluctuations of the 
vacuum, but we don't. The problem enters since, although we do now the explicit form 
of the light quark propagator in the background of a distinct type of vacuum fluctuation 
- instantons (and this propagator has the form we want - poles off the origin), we don't 
know the density of instanton fluctuations in the infrared domain (i.e. where the poles 
are). Although, we know that instantons are not the dominant background fluctuation, 
we are not down on our luck completely, however, we will simply have to develop a model 
of the background fluctuations and fit it phenomenologically from known examples of 
possible duality violations. In the original paper |^], a two parameter model function 
was introduced, and fitted from two sources of duality violations: the semileptonic 
lifetime of the D, and the invariant hadronic mass distribution of r. The model was 
then used to explore possible duality violations elsewhere, e.g. B decays, a s extraction 
from r decays, and more. In the next sections, we review this model. 

10.4 Instantons and the OPE 

Before calculating the actual exponential terms in the real QCD case of D decay, it will 
be useful to review certain universal aspects of the calculation in a toy model. The toy 
model generalizes the actual QCD calculation. Calculating the instanton contribution 
with the toy model, we clearly see that the instanton calculation can give three different 
types of contributions to the transition operator: 

(i) Small-size instantons affect the coefficient functions of the OPE. We are not 
interested in these terms. They do not appear in our calculations below explicitly 
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since our instanton density function excludes them. I will return to this issue when we 
consider the question of the instanton density function in real QCD. 

(ii) Terms proportional to \/{rriQp). They represent the instanton contributions 
to the matrix elements of various finite dimension operators that are present in the 
OPE. In principle, one could calculate each matrix element, pretending that the only 
contribution is from the one-instanton background. Such an attempt leads to results 
grossly violated in nature. 

(iii) The exponential terms. These terms, exp(— m.qp) are what we are after. 
Finding the exponential terms in the case of the heavy quark decay goes just like 

our example of the two-point function above, but with one subtlety - the external lines 
are colored objects, and thus feel the influence of the instanton. This could lead to 
considerable calculational difficulty. Fortunately there is a simplification. To see it, it 
is simplest to consider a toy model where all spins are neglected. The relevant features 
of the toy model translate to the case of real QCD. 

10.5 Instanton Calculation for a Toy Model 

The Lagrangian of the toy model we consider has the form 

L w = hQq<f> + h.c. (10.8) 

which describes the decay of a heavy scalar quark Q into a massless quark q and a 
scalar photon 4>; the coupling h has dimensions of mass. Both quarks are in the spinor 
representation of the color group (SU(2) here). The basic strategy of the instanton 
calculation has been outlined above, here I work out details specific to the heavy quark 
case. 

Consider the transition amplitude: 

T = 2i\b < >= 2^ <Dl I d4xiT i L w(x)Lw(0)]\D >, (10.9) 
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Figure 10.2: The diagram representing heavy quark decay in a toy model in the instan- 
ton background field. 

see figure 10.2. 

For our choice of the toy model Lagrangian, 

f = i J Q(x)S(x,0)Q{0)G^(x 2 )e imcVX d 4 x, (10.10) 

where G§ is the propagator of the scalar photon, and S(x, y) is the propagator of the 
massless scalar quark in the external (instanton) field. The propagator of the massless 
scalar particle in the instanton background is known exactly: 

s(x, y) = . 2( l . 2 (i + p 2 /* 2 )- 1/2 (i + p2{T+ 2 X \ iTy) )(i + P 2 /y 2 r 1/2 , (10-11) 

47r z (x — yy x z y z 
where we have fixed the instanton center at z = 0, and 

T = (f, i); r + = (f, -i); T+Tp = 5 a p + r? a/ 3 7 T 7 , (10.12) 

where r are the Pauli matrices acting in the color subgroup. 

We now take a closer look at the final state quark propagator, rewriting it using the 
Feynman parametrization, 

1 



S(x,y) = - f 1 da[a{l-a)]- l l 2 — 
ir Jo all 



(1 — a)(x — y) 2 + p 2 + z 1 



where 

z = z - xa - (1 - a)y. (10.14) 
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In this form, the analytic structure of the propagator is clearer - we can now easily pick 
up the pole off the origin in our integrations, and avoid singularities at the origin that 
give contributions which we are not interested in. We want the pole: 

x 2 = ±j- -. (10.15) 

a(l - a) 

After doing the integration over time, at m c p S> 1, the remaining integrations are nearly 
Gaussian, and run over narrow intervals, 

& „ A; ( Q _ \f „ _L; {z _ J)2 „ _P_ (10 .16) 

m c ' V 2' m c p y 2' m c y 1 

Thus, one performs the remaining integrations merely by evaluating all of the pre- 
exponential factors at the saddle point. In particular, consider the heavy quark external 
line. The heavy field Q(xo,x) can be written in the leading order as 

Q(x ,x) =Te i fo° A ^ dT Q(O l x) + 0(l/(m cP )) = U(x)Q(0,x) + O(l/(m c p)), (10.17) 

where U (x) is some complicated function (just the exponentiated, integrated instanton 
gauge potential). In principle, the complexity of the function U(x) makes for an ex- 
tremely difficult job of evaluating the transition operator. Fortunately, we are lucky, 
and in the saddle point approximation, U(x) is just equal to l! The heavy quark field 



enters at distances x ~ \J p/m c <C p and, therefore the transition operator is finally 
proportional to Q(0)Q(0). Collecting all remaining factors, 

T(fco) = h 2 Q(0)( G<t>{ 2 ~^ 2) ) J dad 3 ^S e -W(P 2 +- 2 )/("(i-")+- 2 Q(0) (10.18) 
where 

G,(x 2 ) = ^ (10.19) 

is the free scalar propagator. Performing the remaining Gaussian integrations, and 
taking the matrix element (note for scalar particles, < Hq\QQ\Hq >= M^/toq) we 
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finally arrive at 

(2im Q p) 

T = -h 2 — p-, (10.20) 

which gives us a contribution to the total width, 



„t . . 9 sin(2niop) 



Note the oscillatory nature of the result upon continuation to Minkowski space. In ap- 
plications, we set sin(mQ/>) = 1, since we use our model only as an approximate estimate 
of duality violations, and since, anyhow, the suspicion is that the strong dependance of 
sin(niQp) on the heavy quark mass, or position of the fixed-size instanton is somewhat 
artificial. The issue of the the remaining p dependance in the pre-exponential will be 
dealt with in the next section. 



10.6 The exponential contribution for 

We now proceed to the actual calculation of the case of real D decay in the instanton 
model. Let me outline the treatment. First, we write the transition operator in the 
instanton background field. It has the same form as equation (|5.2| ). 

f = gi^LL 2 J Q(x)T^S({x^),z,p)T v Q{Q)l^(x)e^ mvx U A x d A z dp D(p). (10.22) 

The Green function of the light quark, S((x, 0), z, p) is expanded in powers of m q , and 
has the form: 

S(x, y) = -—P (x, y) + G{x, y) + m q A(x, y) + 0{m 2 q ). (10.23) 
m q H 

Technically, difficulties with the instanton model develops at the first step: the 
Green's function of the light quark in the background of one-instanton is not defined 
since the Dirac operator has a zero mode (zero modes are denoted above by Pq, and 
regulated by m q ). However, since the weak amplitude we consider contains left-handed 
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quarks only, the problem of zero modes dissappears completely, albeit somewhat artifi- 
cially. Hopefully, this not fully self-consistent procedure will work satisfactorily enough 
for our purposes. Of course, we have no right, in general, to believe that the one- 
instanton contribution will give us good phenomenology, still, since we scale our model 
from other duality violations, we can hope then that the duality violation hierarchy 
from decay to decay is captured correctly. 

Proceeding with the rest of the calculation see ref. 
We arrive at the result: 

Id = -^r -^rgSin(2m Q p)D( /3 ). (10.24) 

Now, the real question of QCD dynamics enters in full. We need to integrate over the 
instanton size p, and an explicit representation of D{p) is required. The calculation of 
D(p) was first undertaken by t'Hooft |48[j . This instanton density function, however, is 
not what we want. In t'Hooft's calculation, only small-size instantons can be considered. 
His instanton density function took the form: 

D(p)= const. (MQCD) b (10.25) 

where b is the first coefficient in the Gell-Mann Low function. Integrating over this 
instanton density function leaves us not with the exponential quantity we are after, but 
instead power like contributions. This contribution was discussed above. Small size 
instantons are hard fluctuations. Taking them into account should in no way reflect 
the contribution of the series of higher order operators - yielding an exponential term. 
Indeed, considering the small-size instantons we are outside the validity of the standard 
heavy quark expansion (HQE). The standard HQE requires the decomposition of the 
heavy quark field in the form Q{x) = exp^irriQV^x^Qlx) which in the hard instanton 
background becomes inapplicable, as well as the statement that heavy quark spin effects 
are suppressed by 1/rriQ, and so on. 

To get the exponential contribution coming from instantons We need the density 



40 1 for an exhaustive discussion, 
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function in the infrared but, of course, it cannot be calculated there. We thus model 
D(p) with the simplest idea we can think of: a fixed size instanton: 

D(p) = N5(p-p ). (10.26) 

Actually, there is some reason to believe that this density function at least represents 
the character of the the soft instanton fluctuations. Results from the instanton liquid 



model [44 1 show that the density function should have a sharp rise and steep fall-off. 
Thus, we can hope at least that our model correctly captures the essence of the real 
QCD instanton density function. 

10.7 Numerics of the instanton model 

After integrating over the instanton density function, we have the results of our duality 
model for the semileptonic decay expressed in terms of two free parameters: 

T J sl = -To^N 96?r sin(2m c p ) = Fq^N ^ (10.27) 
3 [m c p f 3 (m c/ 9 r 

Let me repeat that since the value of sin(2m c/ og) is sensitive to how close the argument 
is to ?i7r, a very model-dependant feature, we thus set sin(2m c /9o) = 1, and take the 
absolute value of the expresion, taking the conservative point of view that our model 
at best determines an uncertainty due to duality violations. 

Now, we are at a crossroads. Unfortunately we cannot test whether the model we 
have developed is capable of giving us a phenomenologically acceptable result. The de- 
pendance of the result on two unfixed parameters requires us to extract some knowledge 
of duality violations from some other processes. At this stage of developement of both 
theory, and experiment, examples are rare. An obvious place to go hunting is where 
hadronic processes have small energy releases. In fact, there is one place - the invariant 
mass spectrum of hadronic r decays - where we actually do see some preliminary signs 



of duality violations. The issue is considered in detail in ref. [40], here I just sketch the 
results. 
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Figure 10.3: Plot of R( v a \E) with Euclidean and Minkowski truncation errors 
sketched on top of the data. 
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Figure 10.3 is the experimental plot of the quantity R\ V ~ A >(E). From the plot, it can 
be seen that there is, at the tail of the distribution starting at about 2 GeV, evidence 
of some kind of oscillation ~ 10%. This oscillation is exactly the behavior predicted by 
the Euclidean exponential error when continued to the Minkowski domain. In ref. , 
this oscillation was used as an input to determine one of the two free parameters of the 
instanton density function. 

Because other sources of duality violations are scarce, we can fit the remaining 
parameter of the model to the ~ 50% duality violating contribution required in the 
semileptonic D lifetime, and then use the model to predict violations in other processes. 
For example, calculating a smattering of B lifetimes we end up with the result: 

|r 7 (6^ ccs)' 



To(b — ► ccs) 



0.006 (10.28) 



and 



Ar 7 _ AT 1 AT 1 

— — (b — > ccs) ~ 2— — (b — > cud) ~ 5— — (b — > uud) 

Fo r r 

AT 1 AT 1 
~ 16-— (6-> dv) ~ 75-— (&-► s-y) 

Fo r 
Ar 7 

~ 300— - (6 -» ulv), (10.29) 
Fo 

thus, the deviations from duality in B decays are negligible. 

Let me briefly mention one other prediction the instanton model makes with, unlike 
in the case of the B decays, more observable implications. If the approach above is 
applied to the hadronic r width, the deviations from duality are estimated as 

AF 7 (r -» hadrons) ^ Q Qf . 
Fo(t" — > hadrons) 

while this seems like a relatively small result, this 5% uncertainty in the width translates 
into a ~ 30% uncertainty in the value of a s (m T ). This is an interesting result in that 
the low energy determinations of a s currently disagree with measurements at the Z 
peak - except for the case of a s extraction from r decays. Shifman p2| points out that 
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this mismatch between low and high energy a s extraction could be a prompt of new 
physics. The model outlined above should at the very least make those who extract 
a s from r widths a bit nervous, and thus raise an eyebrow at the possibility of new 
physics being seen. At the very least, here the strong case for large duality violations in 
m c where the energy release is ~ 1.40 GeV should make one suspect at least minimal 
violations in r where the energy release is ~ 1.77 GeV. 

Let me make one more closing remark concerning the semileptonic decay b —* ctv. 
This decay was recently measured by [49||. The theoretical prediciton, including 1/m 2 



effects is found in [43]. In the decay, there is an energy release of ~ 1.6 GeV, thus, 
it is tempting to consider that there is perhaps a duality violation on the order 5 — 
10%. Currently, the numerics on both the theoretical and expermental side are not 
trustworthy enough to make an observation of duality violation, but it very well may 
be a place to hunt for duality violations in the future, and thus possibly supply evidence 
for the utility of the instanton model. 



Chapter 11 



Conclusion 

"Physics is rich, but life is richer" - Emil Akhmedov 

We have seen that the physics of the D meson decay is rich indeed! Calculation of the 
total decay width requires all of our theoretical might, and then some. OPE power 
corrections, normal perturbative terms, duality fluctuations - what a problem - and we 
are still evaded. In the end, we found that, at least at present, the semileptonic width of 
the D meson is at best a theoretical kitchen for the study of duality violations. This in 
itself is not a complete dissapointment. Many predictions from other hadronic decays 
concerning CP violation, a s extraction, etc. will come in the future, and almost all, 
due to their Minkowskian nature, will come with the built in assumption of duality. To 
wit, if we want to go from predictions made by QCD without some flip assumption, and 
instead a well investigated one, we need to study D and r decays to better understand 
the assumption of duality. It is not ironic, but rather somewhat expected that in 
our duality investigations, we here too come up against the problem of understanding 
QCD dynamics in the infrared. Unfortunately, it seems that this is a problem which will 
continue to defy solution for some time, and thus, it makes the job of approaching QCD 
physics from the phenomenological side outlined in this thesis all the more important. 
It seems that in QCD, the answer to Fermi's question is still yet to come. 
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